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1. Introduction and statement of results 

We give pole free strips and estimates for resolvents of semiclassical operators which, on 
the level of the classical flow, have normally hyperbolic smooth trapped sets of codimension 
two in phase space. Such trapped sets are structurally stable - see §1.2| - and our motivation 
comes partly from considering the wave equation for slowly rotating Kerr black holes, 
whose trapped photon spheres have precisely that dynamical structure - see ^ From the 
semiclassical point of view an example to keep in mind is given by 

P{z) = -h^A + V{x)-l-z, VeC^{W;R), 

with the classical flow described by Newton's equations: 

The incoming and outgoing tails, T±, and the trapped set, K, are defined by 

r± = {(x,0 : 3M, \ip\x,0\ < M ,t ^Too} , K = T+nT.. 

As explained in ^ it is important to consider more general families of operator pencils. 
The general assumptions will be given in §1.1 but the result is already non-trivial in the 
case presented above: X = M", and P{z) = P — z, P = —h'^A + V{x) — 1. 



Theorem 1. Suppose that P{z) is a family of operators satisfying the assumptions in {1.1 



with a trapped set K which is smooth and normally hyperbolic in the sense of ^1.2 and 
contained in Ui (s X. 

If the symbol of d^P{0) is strictly negative near p-\0) n T^^X and W E C°°{X;R) 
satisfies 

W>0, W\u,= 0, W\x\u,= l, 
where 7i{K) G Ui G U2 X then there exist 6q, z/q > such that for \z\ < 60 we have 

{1/ Iraz , Imz > , 

h-Hogil/h), lmz = 0, 
, Imz > —vqH , 

and in particular, z i— )■ {P{z) — iW)~^ is holomorphic in {\z\ < 6q , Im^; > — z/q/?-}. 
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This result is related to the general principle in scattering theory which in mathematics 
goes back at least to the work of Lax-Phillips and Morawetz: the nature of trapping of rays 
is related to the distance of resonances, which is to say poles of the analytic continuation of 
the resolvent, to the real axis. That in turn is related to energy decay, local smoothing and 
other properties of the propagators. The closeness of these resonances to the real axis is in 
particular related to the stability of the trapped trajectories, with stable trapping giving 
rise to resonances close to the axis — heuristically, these are close to being eigenvalues. 
By contrast, trapped orbits near which the dynamics is hyperbolic leads to resonances 
bounded away from the axis - see for a general introduction. In [29] and [31] a gap 
was established when hyperbolic trapped sets are fractal and a certain topological pressure 
condition is satisfied. In Theorem [T] the trapped set is smooth and has the maximal 
dimension. We assume that the flow is r-normally hyperbolic for every r on this trapped 
manifold in the sense of Hirsch, Pugh, and Shub ^28j and Fenichel [21] . That assumption 



is structurally stable - see §1.2 



The proof of Theorem [T] is based on a positive commutator argument with an escape 
function (4.6) in a slightly exotic symbolic class described in § ^3.2 3.4 A similar logarith- 



mically flattened escape function for more complicated (fractal) trapped sets was used in 
[39] . For the semiclassical analysis near closed hyperbolic orbits similar escape functions 
were used by Christianson in [12] and [13]. In a way, the situation here is simpler as we 
assume that the trapped set has codimension two. However, following our arguments might 
simplify the treatment of closed orbits as well. 

In Theorem [2] in ^ we present a closely related result for resonances. For operators 
P{z) = —h'^A + V{x) — 1 — z with V{x) holomorphic and decaying in a conic neighbourhood 
of MJ^ in C" (in fact, for a larger class of operators with real analytic coefficients in - [23]) 
a more precise resonance free region was obtained by Gerard-Sjostrand [26j. The novelty 
in Theorems [T] and [2] lies in the resolvent bounds and the applicability to C°° coefficients. 
The estimates in microlocally weighted spaces of holomorphic functions in [26j do not 
immediately imply polynomial bounds in h, in the resonance free strips. For more recent 
results involving scattering with hyperbolic trapped sets we refer to [I], [3], [SI], [22], [S3], and 
references given there. 

As examples of immediate applications of Theorem [Tj we give the following corollaries 
which follow immediately from the results of |15j : 



Corollary 1. Suppose that X is a scattering manifold (that is a manifold with an asymptot- 
ically conic metric) and — is the non-negative Laplace- Beltrami operator on X . Suppose 
that the trapped set for the geodesic flow on S*X is normally hyperbolic in the sense of 
{1.2. If r{x) = (1 + (i(a;, xo)^) 2 , where d{x,xo) is the distance function to any fixed point 



xo G X , then for A > 1 



A 
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This implies local smoothing for the Schrodinger equation with a tiny loss of regularity: 

Corollary 2. Under the assumptions of Corollary\l\ we have the following estimate valid 
for any T > (large) and e,6 > (small): 

T 

||r-5-5exp(2tAg)M||^i_^^^^ dt < CT,eAu\\l2(^x) ■ 

Based on this, and assuming that the curvature of the asymptotically conic manifold is 
negative (in every compact set), the results of [7] show that Strichartz estimates hold with 
no loss at all. 

Our motivation for considering this geometric set-up comes from the Kerr black hole. 
This is a family of Lorentzian metrics which solve the Einstein equations and describe 
rotating black holes. We refer to [H] for a survey of mathematical progress on the wave 
equation for these metrics and to [12] for some more recent results and references. In 
the physics literature the decay of waves has been studied in terms of quasinormal modes 
which are the analogues of scattering resonances in this setting - see [30] for a physics 
introduction and [3] for a recent mathematical result which provided an expansion of waves 
in the Schwarzschild -De Sitter background in terms of resonances. 

Obstructions to rapid energy decay occur, heuristically, due to separate mechanisms at 
high and low frequencies. At high frequencies it is expected that the geometry of the 
trapped set plays a key role and it is on this geometry that we focus our attention. As 
recalled below, the trapped set of Kerr is indeed an r-normally hyperbolic manifold (within 
the energy surface) for all r, diffeomorphic to T*S'^ (or S*S'^ if we restrict to fixed energy). 
It is thus of interest to explore the limits placed on exponential local energy decay by this 
trapping mechanism, and this is exactly the role of resonances. That is to say, as the Kerr 
metric is stationary, we may Fourier transform away the "time" variable, and try to study 
the poles of the putative analytic continuation of the resulting stationary operator across 
its continuous spectrum. This motivates considering general operator pencils P{z) in place 
oiP-z. 

In the case of Kerr the principal obstacle, compared to the Schwarzschild analysis [3] is 
the failure of ellipticity of the stationary operator P{z) near the event horizon of the black 
hole, within the so-called "ergo-region." This failure reflects the failure of our timelike 
Killing field (with respect to which we have Fourier transformed) to be timelike in the 
region in question. Thus, we reduce our question to a simpler model problem by cutting 
away the ergo-region. To do this, we modify our stationary operator by considering only 
the form of the operator near its trapped set, and then adding a complex absorbing potential 
to damp waves propagating outward from it. We then consider the complex eigenvalues 
of the resulting non-self- adjoint operator as a proxy for resonances. Such a construction is 
rigorously known to approximate resonances in certain cases [40j. Thus Theorem [l] yields 
a gap in the spectrum of the operator P{z) — iW near the real axis, at high frequency 
(i.e. in the semiclassical limit). Recently, a meromorphic continuation of P{z)~^ and a 
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rigorous definition of quasinormal modes for Kerr-De Sitter black holes have been obtained 
by Dyatlov [IH]. 

Our paper is concerned only with the analysis near the trapped set. Unlike in most 
other mathematical works on Schwarzschild and Kerr black holes - see for instance [8], 
[3], [171 [18j, l22i |23], pi], [13] — this analysis of the trapped set does not use separation 
of variables and properties of the Regge- Wheeler potential. It is carried out in a way 
applicable to the perturbations of the metric. The structure of the trapped set does not 
change under those pertubations but one cannot separate variables anymore - see the end 



of ^1.2 and q2^for more details. 



To indicate how the local results near the trapped set can be used to obtain energy decay 
we present Theorem [3] in [|5j Here is its simplest version: 

Corollary 3. Suppose that X = Xq U (M" \ 5(0, i?)) U ■ ■ ■ U (M" \ 5(0, R)), where Xq is a 
smooth compact Riemannian manifold with boundary, with the metric g equal to the usual 
Euclidean metric in the infinite ends, M" \ B{0,R). If n is odd and the trapped set for 
the geodesic flow on S*X is normally hyperbolic in the sense of ^1.2 , then the local energy 



decays exponentially: for any e > there exists a = a{e) > 0, such that if 

{df - Ag)u = , u\t=o= ui , dtu\t=Q= ui , supp Uj C U d X , 
then for any V X we have 

(1.2) / {\uit,x)\^ + \dMt,x)\^) dx < Ce--*(||«o||W + , 

Jv 

where C depends on U,V , and e. 



Comments on notation. For a set A we denote by neigh (A) a small open neighbourhood 
of A. For V a Banach space, / = Ov{g) means that \\f\\v < C\g\, with the similar notation 
for operators: T = 0{g) : V — )■ W, means ||TM||vy < Unless specified by a 

subscript C denotes a constant the value of which may vary throughout the paper. The 
notation a < b means that a < Cb. 

Acknowledgments . The authors gratefully acknowledge helpful conversations with Kiril 
Datchev, Stephane Nonnenmacher, Clark Robinson, and Amie Wilkinson; Semyon Dyatlov 
and Andras Vasy provided helpful comments and corrections to the manuscript. This work 
was partly supported by NSF grants DMS-0700318 (JW) and DMS-0654436 (MZ). 

1.1. Global assumptions on P{z). We make abstract assumptions on P{z) in order to 
allow very general end structures. The assumptions are in some sense the reversal of the 
black box assumptions of [TU] and [37]: we specify the operator in the compact interaction 
region but allow an almost arbitrary structure outside. That is natural since we are adding 
the complex absorbing potential. Many results about resonances can be rephrased in this 
setting. In some cases they can then be "glued" to obtain global results as was done for 
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scattering manifolds in [15]. Some infinities appear remarkably resilient to that approach, 
in particular the ends of conformally compact, that is asymptotically hyperbolic, manifolds. 
However, we expect that the Kerr metrics can be "glued" to our local construction. 

For a concrete example of operators satisfying the abstract assumptions presented here 
see ^ 

We consider a holomorphic family of operators, 

z^P{z), zeDiO,5i), 

depending implicitly on the semiclassical parameter h. These operators act on Ti, a complex 
Hilbert space with an orthogonal decomposition 

where Xq (e X is an open submanifold of X with a smooth boundary. 

The corresponding orthogonal projections are denoted by IqU and liu respectively, where 
u eT-L. The operators 

p{z) -.n^H 

with the domain V, independent of z (and of the implicit parameter h), and satisfying 

\V = H\Xo), d,P{z):V^n, 
see [37] for a more precise meaning of the first statement. 
We also assume that 

(1.3) \P{z)u = Po{z){u\xo), foTueV, 

where Po{z) G \E'^'°(X), for real values of z, Po{z) is a formally self-adjoint operator on 
L'^{X) given by 

P = p{x,hD) + hp,{x,hD-h), p,{x,hD)e^j:\X), p{x,0>{OVC-C; 



see |:3.1 for the definition of the classes of operators, and for the conditions on X. 

We assume that P{z) is self-adjoint for 2; G M, and that P{z) is holomorphic for 2; G C, 
\z\ < 61. Hence, 

P{z) = P{z)\ \z\<6i, {P{z) ~ i)-^ : n — yV, Im;z = 0, 1^1 < 5i. 
This implies boundedness in a complex neighbourhood, since P{z) — P(Re z) = 0{ \ Im z\) : 

(1.4) {P{z) - t)-^ : n — yV, for|z|<52. 

The assumption that 

1q{P — i)^^ : Ti — > Ti is a compact operator, 

and estimates in ^4.1 imply that {P{z) — iW)^^ : "H — )■ "H is meromorphic in D{0,6i). 
However we do not make this assumption and prove the estimates on the resolvent directly. 



6 J. WUNSCH AND M. ZWORSKI 

As stated in Theorem [T| we further make local assumptions near the trapped set as 
follows: the symbol of d^P{0) is strictly negative near p-^(O) n Tj^^X, and W G C°^(X; M) 
satisfies 

W>0, W\u,= 0, W\x\u,= l, 
where n{K) G Ui G U2 X . Our dynamical assumptions near K follow in the next section. 
Finally we will consider the operator with complex absorbing potential given by 

P{z) - iW 

where we define the operator W by 

\Wu = W{x)u\xo 
with W{x) a smooth function equal to on Ui and 1 on Xq\U2, and 

liWu = liu. 

1.2. Dynamical assumptions. We now discuss the dynamical hypotheses for Theorem[T] 
We first state the minimal hypotheses needed for the proof of the theorem to apply. 

Let if^ denote the flow generated by the Hamilton vector field Hp. Let r denote the 
distance function to a fixed point in X and locally define the backward/forward trapped 
sets by: 

r± = {p G 7r-\U2) : lim r(^*(p)) ^ 00}. 

Let = F-i- np^^(A). We can then define the trapped set 

K = T+r\T_ 

and let Kx = Knp-^{\). 
Dynamical Hypotheses. 

(1) There exists 6 > such that dp ^ on p~^{X) for |A| < S. 

(2) F-i- are codimension-one smooth manifolds intersecting transversely at K. (It is not 
difficult to verify that T± must then he coisotropic and K symplectic.) 

(3) The flow is hyperbolic in the normal directions to K within the energy surface: there 
exist suhhundle^E^ of Tx^{T^) such that 

TKji=TKx®E^, 

where 

and there exists 6* > such that for all\\\ < 5, 
(1.5) \W{.v)\\ < Ce"^l*l||t;|| for all v e E^ , ±t > 0. 

^The bundles may of course depend on A but we omit this dependence from the notation. 
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These assumptions can be verified directly for tlie trapped set of a slowly rotating Kerr 
black hole (i.e. when a is small) but they are not stable under perturbations, hence do not 
obviously apply to perturbations of Kerr. However, we will show that Kerr in fact satisfies 
a more stringent (and well-studied) hypothesis that is stable under perturbation, and that 
implies the Dynamical Hypotheses above. In particular, the standard dynamical notion of 
r -normal hyperbolicity implies items ^ and (|3]), and is stable under perturbations, modulo 
possible loss of derivatives: 

Recall that the flow in the energy surface p~^{X) near Kx is eventually absolutely r- 
normally hyperbolic for every r in the sense of [281 Definition 4] if its time-one flow is a C" 
map preserving a manifold JCx (which a priori need only lie but is then automatically in 
C") such that for all p G Kx, there exists a splitting of the tangent bundle into subbundles 
stable under the flow 

T,p'\X) = T,Kx © E; © E; , dv'^{E^) = Ej(^) , 

and for each r G N there exist > and C > (both depending on r) such that for t > 0, 

(1.6) 



sup \\dipl\TKj <Ce-'''' inf \\dip-^'\E+f' 



inf \\d^~'\TK,\r >C-^e'''^ sup W,\e-\\ 

with ||«|| some (indeed, any) fixed Finsler metric. This assumption thus entails not merely 
that there is expansion and contraction in the normal direction to K but also that this ex- 
pansion/contraction is considerably stronger than any expansion and contraction occuring 



in the flow on K itself. We remark that one may easily check that (1.6) is stronger than 
( 1.5[ ) by noting that since y)* are all diffeomorphisms, fixing a Riemannian metric gives 



sup ||(i(y9*(p) II > 1 



for all t, hence, for instance the first line of (1.6) gives the estimate (1.5) for the bundle 

E+. 

We may replace hypotheses ^ and ^ with the assumption that for |A| < 5, the trapped 
set Kx has the property that the flow near it in p~^{X) is eventually absolutely r-normally 
hyperbolic for every r. The existence of manifolds T± tangent to E^ and satisfying the 
Dynamical Hypotheses, as well as the structural stability of these assumptions, are clas- 
sical theorems of Fenichel |2T] and Hirsch-Pugh-Shub [28]. The resulting perturbed sta- 
ble/unstable and trapped manifolds are only finitely differentiable in general, as r-normal 
hyperbolicity for each r is the structurally stable property, and this only entails regu- 
larity; on the other hand this r can be chosen as large as desired. While we stated the 
theorems above with C°° hypotheses for simplicity, it is manifest from the proofs that the 
hypotheses could be reduced to insisting that T± be in for sufficiently large K, hence 
those results apply to the perturbed trapped sets arising here. 
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Thus once we show in the following section that the trapped set for Kerr satisfies the 
r-normal hypcrbolicity assumptions, we will know that perturbations of Kerr continue to 
satisfy the Dynamical Hypotheses, with as much differentiability as is required. 

2. Trapping for Kerr black holes 

The hypotheses in the preceding sections are motivated by the example of the slowly 
rotating Kerr black hole. In this family of examples, describing the geometry of a rotating 
black hole, the structure of the trapped set is as described above, while the global structure 
of the spacetime is more complex. The proof that the Kerr trapped set is r-normally 
hyperbolic might be a new contribution. 

We now recall the Kerr geometry, and verify that the hypotheses from the preceding 
section hold in a spatial neighbourhood of the trapped set, at least for small values of the 
parameter a describing the angular momentum per unit mass of the black hole. 

The Kerr metric is a metric given in "Boyer-Lindquist" coordinates by 
with 

= + cos^ 9, 

A = - 2Mr + a^. 
We study this metric on R x (r+, oo) x S'^ with 

in this region, outside the "event horizon" r = r+, the metric is a nonsingular Lorentzian 
metric. The parameter a G [0, M) is the rotational parameter (angular momentum per unit 
mass), and M is the mass. When a = we have spherical symmetry, and the Kerr metric 
reduces to the Schwarzschild metric. 

The d'Alembertian in the Kerr metric is given by 
□ = ( a sm e)dt + -^dtd^ + ( A ~ ^^^^ ~ ~ 

Thus, setting Dm = 0, if is of the form e*^*V£;(r, 6, <^), we find that ve satisfies PeVe — 0, 
where Pe is given by 

-E\^ ^^-a^sin^e)+iE——d^-{-- + ^^)dl-d,Adr-—desmede. 

^ A ' A ^ A sm 0' ^ smU 
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Setting E = {1 + hw)/h (and dropping the subscript on v) we have 

(2.1) (^{l + 2hw){~^^-^^ + a^sm^e)-{l + hw)^^hD^ 

+ ( - ^ + -A7^)i^D^f + ihDr)A{hDr) + -^{hDg) sme{hDe) + 0{h^)]v = 0. 
Thus, if we set 

(2.2) P={- ^ ^ ' +a^smH) - ^^{hD^) 

+ ( - X + -At,) ihD^f + (m,)A(m,) + —AhDe) smOikDe) + 0{h^) 
A sm a sm a 

and 

Q = 2(^ - sin^ 9) + -^{hD^), 

we are deahng with the equation 

(P - hwQ)u = 0. 

The operator P has disagreeable asymptotics near the ends r = r+, oo, however; we thus 
choose to multiply our equation through by A/r^. Thus, we let P = (A/r^)P and Q = 
(A/r^)Q, so that 

(2.3) P^(-<^ + ^s,„^«)-^^(.Z,.) 

+ ( - ^ + ^i^X'-.O^)' + ^(ftAOACftA.) + ^g(I^Oe) sin9(AD„) + O(ft^) 

and 

Q = 2(^^:^ - sin^ e) + —rihD,), 

and we are now interested to solutions of P{z)u = 0, with 

P(z) =P-zQ, z = hw. 

We are in the situation covered by Theorem [T] provided that we can verify the hypotheses 
on P and P'{0) = —Q. We note that P and Q are now self-adjoint with respect to the 
volume form 

4 

r r. — r 

— y'\g\drd6 dip. 



To see this, we write 

^(_iA^ 
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where P' is our original, formally self-adjoint operator □, applied to functions independent 
of the t variable (i.e. on the quotient of the spacetime by the dt flow); the terms are 
self-adjoint by axial symmetry of g. 

The hypotheses are, we claim, satisfied in a subset {r > tq} (for some tq > r+) that 
includes the trapped set and the r — )■ +00 end. The hypotheses are not globally satisfied, 
however, owing to the structure of P near the event horizon: not only is this end not asymp- 
totically Euclidean, but the operator P is not even elliptic in a uniform neighbourhood of 
r = r+ : inside the "ergosphere" where 

1 

P is not elliptic (i.e. the Killing vector field dt for the Kerr metric fails to be timelike). 
Thus we do not at this time know how to fit the global structure of the Kerr metric into 



the assumptions made in [1.1 for the moment we would instead have to consider a Kerr 
metric glued to a Euclidean end in place of the r — )■ r+ end. 

In what follows, we verify that the structure of the Kerr trapped set, at least, is of the 
desired form. Letting 

^dr + ad6 + P dip 

denote the canonical one-form on T*X, we find that the semiclassical principal symbol of 
P = (rVA)P i^ 

(2.4) P = + " + (-2^ - A )^ - ( A " ^ ^) 

and the Hamilton vector field is given by 



{a{aP + 2Mr) - /3 A csc"^ 9) 

aT 

+ {13^ cot 6 csc^ 6 — 0? sin 6 cos 6) da 



(2.5) {1/2)H = ^Adr + adg - ^-d^ 



+ ((M-r)e + jd^. 

We note (following Carter pj) that the quantities 

p, j3, and JXf = + (asin6' 

^ sm 

are all conserved under the if -flow, and in involution, both on and off' the energy surface 
{p = 0}. 



4ii our analysis of the null bicharacteristics, we study the operator (r'^/A)P, which of course has no 
effect on the dynamics on K\. 
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Under the H-fiow, for each fixed (3, the sets of variables {6, a) and (r, ^) evolve au- 
tonomously, with describing a conserved quantity in the {6, a) plane. This demonstrates 
that the motion in the {6, a) variables is periodic. Also, 

is conserved and (for /3 fixed) dependent solely on (r, ^). This last observation means that 
in fact under the rescaled fiow, generated by {1/2 A)H, the quantity 

_.2 . a2/32 + 4Mar/3 + (r2 + a2)2 
is constant. For a = 0, this quantity is simply 



A' 



4 

The "potential" — ^ has a nondegenerate local maximum at r = 3M; this is its only critical 
point outside the event horizon. Thus this rescaled fiow tends to r = +oo or r = r+ except 
when r = 3M, where it has an (unstable) invariant set (r = 3M, ^ = 0). More generally, 
for a small, the structure is more or less the same: for each given /3, there is a unique local 
maximum of the potential 

, , , ^ a2/32 + 4Mar/3 + (r^ + a^f 
vi3{r) = 2a(3 

outside r = r+. Thus, the trapped set K consists of a family of orbits on which r = r(/3), ^ = 
0, with r{P) given by the critical point of vp in the exterior of the black hole. The invariance 
of p and (3 on the four dimensional trapped set r = r(/3), ^ = with coordinates (6*, ip, a, (3) 
yields the desired integrability. (Note that p and (3 are manifestly in involution.) 

To verify the hypothesis (1.5), we note that since the center manifold is given by r = 
r(/3),^ = 0, we need only verify that the fiow in r, ^ is hyperbolic near these points. The 
linearization of this fiow is simply 

A(r; 
B'ir) ' 



where, by (2.5), 

ia(3{M - r) + rA + M{a'^ - r^)) {a(3 + (a^ + r^)) 
= 

The positivity of B'{r) at r = r(/3) is equivalent to the positivity of v4'(r), where 

A{r) = {a(3{M - r) + rA + M{a^ - r^)) {a(3 + (a^ + r^)). 

When a = 0, strict positivity is easily verified at r = r{(3) = 3M; again by perturbation, it 
persists for small a. 
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We note that in the special case of the Schwarzschild metric (a = 0) we can simply 
compute from ( 2.5[ ) that at the trapped set r = 3M, ^ = : 



where primes denote derivatives under the flow generated by {1/2)H. Thus the unstable 
Liapunov exponent under the H-How is Qy/SM. 

For any given /3, let 7^ denote the subsets of M^^^ given by the stable and unstable 
manifolds of the fixed point (r = r(/3),^ = 0). As /3 is conserved under the flow, the 
fibration 

{(r, e, e, ^, «, (3) : (r, ^) e ^f} ^ (r = r(/3), ^ = 0, 9, ^, a, /3) 

gives smooth fibrations of the stable and unstable manifolds of the flow. (The fibration is 
conserved under the flow since 7^ and /3 are.) 

To check the hypotheses on Q = —P'{0), we note that 

cr{Q) + p = ( 7 sin^ ^) + ( X" "I 9~)f^'^ + nonnegative terms. 

A ^ ^ A sin 6 

The first term on the right is bounded below by 

+ aV^ + 2MaV 



A 

while the second is bounded below by 
(2.6) 



A 

hence we obtain the positivity of (j{Q) (hence negativity of -P'(O)) in a spatial neighbour- 
hood of the trapped set, provided a is not too large; recall that for a = 0, the trapped set 
lies over r = 3M, where the latter term in (2.6) is safely positive. 

We now show that the hypotheses of Theorem [T] are indeed satisfied near the trapped set 
not just for the slowly rotating Kerr metric itself, but for smooth perturbations of such Kerr 
metrics. The crucial observation is that for a small, the Kerr metric is r-normally hyperbolic 
for every r, and that these properties are structurally stable, so that an invariant manifold 
diffeomorphic to S*{S'^) persists, with the flow near it remaining normally hyperbolic. We 
recall that the perturbed trapped set may cease to be infinitely differentiable: for any r, 
a sufficiently small perturbation gives a trapped set in but the required perturbation 
size may shrink as r — 00. In practice this need not concern us, as the proof of Theorem [T] 
only uses a finite (albeit unspecified) number of derivatives. 

Proposition 2.1. For a sufficiently small, there exists a neighbourhood of K, such that the 
flow generated by H is r-normally hyperbolic for each r, i.e. satisfies (1.6). Hence, by the 
results of [28], for each r, any sufficiently small perturbation of the Kerr metric also gives 
rise to an r-normally hyperbolic trapped set (in satisfying the hypotheses of [1.2. 
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Proof. We have verified above that dip^\E± satisfies 

\\dipl{v)\\ < Ce-^l*l||t;|| for all v G Ej, ±t > 0, 



for some ^ > 0. To further verify (1.6) we also require estimates on dip \tk ■ Recall that the 
fiow on K is integrable for the simple reason that p and (3 are both conserved (i.e. we only 
use axial symmetry here, not preservation of as well). Fixing the values of p,P foliates 
K into invariant tori on which the fiow is necessarily quasi-periodic. As a consequence of 
the quasi-periodicity, away from any possible degenerate tori, we have action-angle variable 
(Ji, . . . , 4) G M", (^1, . . . G (^^)" such that H = J2 ^jWe,, hence 

d{^\pW{p')f ~ - r^f + [e, - er + - u,{i'))tf < rf(p,p')'(i + {tf). 

Thus, 

\\dv'\TK\\ < C{t). 

Near degenerate invariant tori, this argument breaks down, and could in principle fail 
(e.g. there can be hyperbolic closed orbits on surfaces of rotation). However we claim that 
the same estimate in fact holds globally on K; it thus remains to check it near degenerate 



tori. Restricting p given by (2.4) to the trapped set, where ^ = and r = r(/3), we find that 
dp and dp are linearly dependent only at a = 0, 6 = 7c/2, i.e. at the equatorial orbits. (A 
separate computation shows that orbits passing through the poles, i.e. with /3 = are not 
degenerate, even though the coordinate system employed here is not valid near the poles.) 
Put another way, the functions f3 restricted to the set KDp'^^X) has its only critical points 
along the set a = 0, 6 = tt/2,lp G S^. In the case of the Schwarzschild metric (a = 0), 
there are two values of /3 at which this can occur, ±{E + r^/A) and they are respectively 
maxima and minima nondegenerate in the sense of Morse-Bott. In particular, we may use 
coordinates a, 9,ip on K np~^(A), and for the Schwarzschild case, K = {r = 3M} and 

4 

/3 = ±sin^(A + ^ - a')'^', 

hence at the critical manifold 6 = 7i/2,a = we compute 

/3L = T{\ + 27M2)-i/2, = T(A + 27M'y/', p% = 0. 

This establishes nondegeneracy, which extends by continuity of second partial derivatives 
for the Kerr case when a is small. 

The behavior of an invariant torus in a three-dimensional energy surface near a Morse- 
Bott maximum or minimum of a conserved quantity is well understood (see, e.g. [2J): it 
must be an invariant circle surrounded by nondegenerate invariant tori shrinking down to 
it; in particular, if /3 takes on a maximum values Pm-, along an equatorial orbit, then any 
sufficiently nearby orbit is constrained to lie for all time in /3~^((/3m — e,/3Af)), and this 
is a solid torus x B"^ in the energy space surrounding the equatorial orbit S"^, whose 
diameter can be made as small as desired by shrinking e — > 0. Taking a cross section of 
this solid torus, we observe that the Poincare return map is thus a twist map preserving 
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the value of /3, under whose iterations the distances between points grows hnearly in time. 
Additionally, we of course have = along the flow, so the difference between /3 values 
can grow at worst linearly along the orbit. Thus, we again obtain linear growth of distances 



along the orbit, hence dip \tk grows at most linearly. This implies (1.6) for every r. □ 



We have thus established the dynamical hypotheses for the Hamilton vector field H, 
associated to p. As 2; G M varies, this is not all of the real part of the symbol of P — zQ; 
by structural stability, however, the hypotheses persist for the principal symbol of P — zQ 
for 2 G M sufficiently small. 

Finally, we observe that in the the end of the manifold r — )■ +00, the assumptions on P{z) 
can be routinely verified by use of the semi-classical scattering calculus of pseudodifferential 
operators [H], as P{z) — iW is elliptic in that setting. 

3. Analytic preliminaries 

In this section we recall facts from semiclassical analysis referring to [16] and [20] for 
background material. 

3.1. Semiclassical calculus. Because of our assumptions, except in ^ we will only use 
semiclassical calculus on a compact manifold. Thus, let X be a C°° manifold which agrees 
with M" outside a compact set, or more generally has finitely many ends diffeomorphic to 
W : 

(3.1) X = Xo U Xi U ■ ■ ■ U Xtv , where Xj = \ P(0, R) for j > 0, and Xq m X. 
We introduce the class of semiclassical symbols on X (see for instance [201 §9-7]): 

S''n,k^r^*X) = {ae C°°(T*X X (0, 1]) : {d^d^aix,^; h)\ < Ca,ph'^ {0""~^^^} , 

where outside Xq we take the usual M" coordinates in this definition. The corresponding 
class of pseudodifferential operators is denoted by \E'™''^(X), and we have the quantization 
and symbol maps: 

ah ■■ ^"■''(X) — > S"''''{T*X)/S'"'-^'''-\T*X) , 
with both maps surjective, and the usual properties 

ahiAoB) = aHiA)ahiB), 

(3 2) 

^ ^m-i,/=-i^^) ^ ^ S"'^\T*X)/S"'-^^^-\T*X) , 

a short exact sequence, and 

(Th o Op)^ : 5'"'*^(T*X) ^ S"''\T*X)/S"'~^'''-\T*X) , 
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the natural projection map. The class of operators and the quantization map are defined 
locally using the definition on M": 

(3.3) Op-Mu{x) = J I' a e^^-y'^^/'u{y)dyd^ . 

We remark only that when we consider the operators acting on half-densities we can define 
the symbol map, ah, onto 

We keep this in mind but for notational simplicity we suppress the half-density notation. 

For future reference, and to illustrate the uses of the calculus, we present the following 
application: 

Proposition 3.1. Suppose P G \E'^'°(X) satisfies P = p{x, hD) + hpi{x, hD; h), pi G ^1/^'°, 

p{x,o>{OVc-c. 

(%) Let tfj, G C^{T*X; [0, 1]), j = 1,2, satisfy 

i,^ = 1 m p-\[-j5,j6])) = 1 , supp^, C p-\[-{j + 1/2)5, (j + 1/2)5]) . 
Then there exists Ei G \1/;^^'°(X), such that 

E,oP = I + R,, i?iG<°(X), 

and 

(1 - ^^{x, hD))R, G M/,°^'-°^(X) , hD)E^ G M/,°^'-°^(X) . 

(a) Suppose f G C^(X, [0, 1]), satisfies f = 1 on U C X, U open. Then there exists 
E2 G *^^'°(X), such that 

E20{P~tf)=I + R2, i?2G<'°(X), 

and 

Xi?2, i?2Xe^r'""W, foranyxeC^iX), snppx<^U. 

3.2. Si spaces with two parameters. As in [521 §3.3] we define the following symbol 
class: 

/~\§(H + I/3|) 

(3.4) a G ST'^^'^T^W") ^ \d^dlaix,0\ < Co.ph-'^h^^ I ^ I (O'''"' , 

where in the notation we suppress the dependence of a on /i and h. When working on M" 
or in fixed local coordinates we will use a simpler class 

(3.5) a G 5i (T*M") ^ |5"a| < C„,jv {h/h^ (^"^ ■ 



rm+m' ,m+m' -M,k+k' -M 
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Then standard results (see pQ^j §9.3]) show that if a G ST'"^''' and b G ST '"^ ''^ then 

2 2 

a(x, /i-Dx) ° ^(a;, /i-Dx) = c(x, hD^) with c G ST'^"^ ,rn+m ,k+k ^ 

2 

The presence of the additional parameter h allows us to conclude that 
c = V ^dfaD^b mod ST 

\a\<M 

that is, we have a symbolic expansion in powers of h. We denote our class of operators by 
^?'™'''(T*M"), or in the case of symbols in 5i, 5i . 

A standard rescaling shows that this class of pseudodifferential operators is essentially 
equivalent to the calculus with a new Planck constant h\ put 

(3.6) (x,0 = (V/i)5(x,e), 

and define the following unitary operator on L^(]R"): 

Uh,h^{^) = ih/h)^u{{h/h)-^x). 

The one easily checks that 

a{x,hD^) = U-^^a^jX^,hD^)Uf^ -^, %,^(x,0 = a{{h/h)^x,i)) . 



Clearly a satisfies (3.5) if and only if a^^ G S{T*W'), with estimates uniform with respect 
to h and h. 

We recall [391 Lemma 3.6] which provides explicit error estimates on remainders. 
Lemma 3.2. Suppose that a,b E Si, and that d" = oh^ . Then 

(3.7) c{x,0 = ^—{—a{D^,D^]Dy,Dri)\ a(x,0^(l/,^)fx=y,c=r, +eiv(a;,^) , 

where for some M 

(3-8) X V sup sup (/ii%,5;,,,))^(za(Z})/2)^+i9-^a(x,e)9"^&(l/,^) 

where ^{D) = a{Dx, D^; Dy, D^) . 
As a particular consequence we notice that if a G 5*1 (T*]R"') and b G S{T*W^) then 



k\ 

k=0 
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3.3. The \E'i calculus on a manifold. On a manifold of the type defined in the beginning 



of ^3.1 



we consider the following class Si: 

2 

Si = Si{T*X) = {a e C^{T*X) : dZ,.a = {h/hy^''^^^O{{0~°")} 



where outside of a compact set we use Euclidean coordinates, determined by the infinite 
ends of X. 

We first observe that this class is invariant under symplectic lifting of diffeomorphisms 
of X, constant outside of a compact set. To define \E'i(X) we need to check invariance of 

5*1 (T*]R") under local changes of coordinates. Towards that we have the following lemma: 

Lemma 3.3. Suppose that a G 5'i(T*]R"), Uj C M", j = 1,2 are open, and f : Ui ^ U2 is 
a diffeomorphism. Let x G C^{Ui). Then A2 == x^^^i^^ hD)x = a.x(^' ^^)> where a^&Si, 
«x = XO'X + (h^h^)- For Ai ^= {f^^)*Af*, we have 

Ai = a'"f{x,hD), a/G5i(T*M"), 

and 

(3.10) af{x,i) = X{r\x))a{r\x),'f'{x)0x{r\x)) + 0~,^ {hW^) . 

2 

Remark. It seems important that we use the Weyl quantization. In the case of the right 
quantization 

a\x,hD)u = j^^ j j a\x,i)e'^^-y^^^'^u{y)dydi, 
we have the exact formula 

see [271 (18.1.28)]. The asymptotic expansion [27, (18.1.30)], 

«;(/(^)'^)~ E ^(5^i)(a^,*/'(a:)r7)(/i/^,)V<''-(^)'''>/nx=,, 

px{y) =^ f{y) - f{x)- f'{x){y - x) , 

is valid in our case as an expansion in h only. In fact, due to the second order of vanishing 
of at X, 

and 

id^a^){x,'f'ix)7]) = 0{{h/h)-\''\l\r])-^ . 
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Hence the terms in the expansion are in 

2 

(the term with |a| = 1 vanishes). 

The Weyl quantization will also be important in local arguments in §4.2[ Finally we 
remark that for this class of symbols the improvement in the error occurs only in h when 
the action of half-densities is considered - see [3Sl Appendix] or [2Qi Theorem 9.12]. 



Proof. The statement about follows from Lemma 3.2 For the change of variables we 
consider the Schwartz kernels of A2 = a'^ix, hD) and Ai = aj{x, hD) as densities: 

(3.11) K,ix,y)\dy\ J b e^^^-y'^^/'d^\dy\ , 
which means we seek af such that 

(3.12) Ka^{x,y)\dy\ = Ka,{x,y)\dy\, x = f{x) , y = f{y) . 
We rewrite the right-hand side as by changing variables 



Writing, 



fix) - f{y) = Fix, y){x - y) , F{x, y) = f ( ^ ) + 0{{x - yf) , 
(3.13) ^ ^ 

m + f{y) = fl^^^j+0{{x-yr). 

we apply the "Kuranishi trick" by changing variables in the integral, = F{x,yyC,'- 



{2TThy 



(27r/i)" 



xe^'^^-y'^>/'^d^\F{x,yynf'{y)\\dy\ 

/("^(/(^)^(/'(^)TM-«.,('--^(^-)= 



xe'^^-y'^^/'^d^\f'{{x + y)/2)nny)\\dy\. 
We now observe that 

\n{x + y)/2)\ = \ny)\ + 0{\x-y\) 
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and consequently Kaf{x,y)\dy\ 
1 




The terms 



(hh-Hx-yf) 

~ 3 1 

contribute terms Og (h^h^) to the symbol: we use integration by parts based on 

2 

(x - y) exp((x - y,0/h) = hD^exp{{x - y,C)/h) . 
Similarly, smooth terms of the form Og {\x — y\) give contributions of the form (h^h^). 

Here in dealing with the "big-Oh" terms we use the fact that for h = b{x,y,^) G Si (with 
the definition modified to include derivatives with respect to y), 



{27Thy 



h{x,y,i)e''^^-y'^^l^di 



{2'KhY 



X + y 



where 



b^{x,0=b{x,x,0 + Os^ (h) 



which follows from the standard pseudodifferential calculus and the rescaling (3.6). 
This shows that Kaf{x,y)\dy\ = 



{2'KhY 



(^)'(/'(^)X'^)^^('-"'))) 



J{x-y. 



hence aj can be chosen in the form (3.10) so that this matches Ka{x,y)\dy\. 



□ 



We need one more lemma which shows that away from the diagonal the symbol contri- 
bution is negligible in h (rather than merely in the h sense). This does not contradict the 
rescaling (3.6) which eliminates h, as the distance to the diagonal then grows proportionally 
to h-^/^ (see [201 Theorem 4.18]). 

Lemma 3.4. Suppose that Xj G C^(M"') are independent of h, and suppxi nsuppx2 = 0- 
Ifae 5i(T*M") then 

Xia}"{x,hD)x2 = Os'^s{hn- 



Proof. We can apply Lemma 



3.2 



as in the composition formula for a E Si and h E S 

2 

presented in (3.9): in the composition Xi(^^X2 all terms in the expansion vanish and the 
error becomes arbitrarily smoothing and bounded by /i^, for any N . □ 
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Using Lemmas 3.3 and 3.4 we obtain an invariantly defined symbol map for tlie class 
(X) defined using local coordinates, as in [271 §18-2] (see [20l §E.2] for the semiclassical 
case). The symbol map occurs in the following short exact sequence: 

— y /lUi^i(X) — ^ ^i(X) ^ Si{T*X)/h'2h'2Si{T*X) — ^ 0. 



2 



This means that if we start with a G h^"^Si{T*X) then the operator a^{x,hD) G 
h'^"^ i{X) is well defined and its symbol is determined in any local coordinates up to 
terms in h^^^'ihhSi. We will be particularly interested in the case 

(3.14) a e Sl{T*X) = Pi h-'^S.{T*X) , 

in which case the local symbols will be determined up to terms of size h^h^S]^. 

2 

3.4. Exponentiation and quantization. As in [3S] and [12] it will be important to 
consider operators expG^{x, hD), where G ^ S^. To understand conjugated operators, 

2 

exp(-G""(a;, /iD))Pexp(G"'(x, hD)) , 

we will use a special case of a result of Bony and Chemin [5, Theoreme 6.4] - see [39l 
Appendix] or [201 §9.6]. Because of the invariance properties established in j j3.3| we discuss 
only the case of M" in the next two subsections. 

Let m(x,^) be an order function in the sense of |16] : 

(3.15) rn{x,^) < Cm{y,r]){{x - y, ^ - rj))'^ . 
The class of symbols, S{m), corresponding to m is defined as 

a e S{m) ^ |9^9g^a(x,0l < Capm{x,0 ■ 

If mi and m2 are order functions in the sense of (3.15), and aj G S{mj) then (we put h = 1 
here), 

a^(x, D)a^{x, D) = }f{x, D) , be ^(mims) , 
with b given by the usual formula, 

6(x,0 = ai # 02(3;, 

(3.16) 

= exp{ia{D^i,D^i;D^2,D^2)/2)ai{x^,^^)a2ix^,^'^)\x^=x2=x,e=e=i ■ 
A special case of [5, Theoreme 6.4] (see [39l Appendix]) gives 



Proposition 3.5. Let m be an order function in the sense of (3.15) and suppose that 
G E C°^(r*M";M) satisfies 

(3.17) G(a;,e)-logm(x,e) = 0(l), d^d^^Gix, = 0(1) , |a| + |/3|>l. 
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Then 

(3.18) exp(tG"'(x, D)) = B^{x, D) , Bt e S{m') . 

Here exp{tG^{x, D)) is constructed by solving dtu = G^{x,D)u, u E S. The estimates on 



Bt G S{m^) depend only on the constants in (3.17) and in (3.15). In particular they are 
independent of the support ofG. 

Since m* is the order function exp(t log m(a;, ^)), we can say that on the level of order 
functions "quantization commutes with exponentiation" . 

3.5. Conjugation by exponential weights. Let m be an order function for the Si class: 

/ / ^ ^ 

m(p) < Cm(p') ( — — 

\{h/h)2 

for some A^. We will consider order functions satisfying 
(3.19) mESi(m), — eSi ( — 

This is equivalent to m(p) = expG'(p) with 



'g(p) P-P' 



N 



(3.20) ^^^^^ < C ( ) , = 0{{h/h)-\'^\/' , |a| > 1 

expG(p') - \{h/h)2 ' ' ' - 



Using the rescaling (3.6) we see that Proposition 3.5 implies that 



exY>{sG'"{x,hD))=Ff{x,hD), F,eSi{m'), seR, 

(3.21) 

A = Op^a) , a G Slim') ^ A = exp(sG"'(x, hD))Op'^{ao) , oq G . 
For P G \I^°''^(X) we consider 

(3.22) P,G = (,-'^G^{^^D)/hp^G^{x,hD) ^ g-ad,G»(a=,hD) p ^ $i 

2 ' 

where used Proposition |3. 5 as described above. In particular we have an expansion 



(3.23) ^-sG-i.,kDyHp^sG-i.,HD)/H -^T ^^G-i.,,D)Y P 



i=0 



where 



(3.24) ad^»(,,,^)PG/i/i^-ivl>, 
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3.6. Escape function away from the trapped set. Here we recall the escape function 
from [25; Appendix]. Suppose that U, V are open neighbourhoods of of fl p^^{[—6, 6]), 

fj (£ r (E T*X . 

There exists d G C°°(T*X), such that 

(3.25) Gi\u=0, HpGi>0, -ffpGitp-i([-25,2<5])< C* , -f^pG'ifp-i([-5,5])\y> 1 • 

Since HpGi > 0, Gi is an escape function in the sense of [23]. It is strictly increasing along 
the flow of Hp on p^^{[—6, 6]), away from the trapped set K. Moreover HpG is bounded in 
a neighbourhood of p~^{[—26, 26]). Such an escape function Gi is necessarily of unbounded 
support. 

4. Proof of Theorem [T] 



In ^4.1 4.3 we identify P{z) with Po{z) and assume that u is supported in Xq. In ^4.4 



we will show how the assumptions on P{z) in ^1.1 give a global estimate on the inverse. 
Since we have not assumed that \{P{z) is a compact operator we do noi prove that 

{P{z) — iW)~^ is a meromorphic family of operators. We prove that the inverse exists for 
Im^; > —v^h by direct estimates. 



4.1. Estimates for Im^; > 0. To obtain the first estimate in (1.1) we adapt the proof of 
[38| Lemma 6.1] to our setting. 

For that let ip"^ = ^"'(x, hD), ip e C^{T*X, [0, 1]), be a microlocal cut-off to a a small 
neighbourhood of p~^(0) fl T^^X, and suppose that 

V = {P{z)-iW)u. 

Semi-classical elliptic regularity gives 

(4.1) 11(1 - ^'")u\\ < G\\v\\ + 0{h°^)\\u\\ 



(see part (i) of Proposition 3.1). The assumption that dzP{0) has a negative symbol on 
the characteristic set of p, in the region where < < 1 implies that 

P{z) = P(Re z) -ilmz Q{z) , 

where P(Rez) is self-adjoint and a{Q{z)) > 1/G > near p~^{0) nT^^X. This shows that 

-lm{{P{z) - iW)ij'"u,ij'^u) = lmz Re{Q{z)ilj'"u,'^'"u) + {Wip'^u^ip'^u) 

^^■^^ >Imz {W^P'^uf/C - 0{h'^)\\uf) + {W^lJ'"u,^|J'"u) , 

where we used the semi-classical Garding inequality (see fl6[ Theorem 7.12] or [20l Theorem 
4.21]). We also write 

lm{P{z)u,u) -lm{P{z)'^'"u,'^'"u) =lmz {{Q{z)u,u) - {Q{z)ip'"u,ilj'"u)) 

= Imz 0(1)||(1 -^'")u]|]|m|| 

= lmzO{l) {\\v\\\\u\\ + Oih'^^uW^) , 
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where we used elhptic regularity (4.1) in the last estimate. Then, applying (4.2), 
llMllllfll > — lm{{P{z) — iW)u, u) 

= -lui{{P{z) ~ iW)ij'"u,ij'"u) -Imz Oil) {\\v\\\\u\\ + 0{N^)\\u\\'^) 

+ {{W - ^Ij'"W^Ij'")u,u) - 0{h)\\u\\'^ 
>\mz {Wij'^uf /C - 0{l)\\v\\\\u\\ -0(/i)||nf) . 

Here W — ip'^Wip^ > —0{h) follows from the semi-classical sharp Garding inequality. 

For small Im^; the term ||f || ||-u|| on the left hand side can be absorbed in the right hand 
side, and by adding Im2;||(l — ip^)u\\'^ to both sides we obtain 

Im2;||M|p/C < + 0{h) lmz\\u\\'^ , 

and that gives 

„ „ C .. „ 
\ u\ < 



Im^; 



Combined with the estimates in ^4.4 this proves 



\\{P{z)-tW)-'\\L2^L'^<-^, for lmz>0, \z\ < 6o . 

Im z 

4.2. Estimates on the real axis. In this section we will use a commutator argument to 
obtain an estimate on the real axis. In fact, this bound automatically gives holomorphy of 
{P{z) - iW)-^ in Im;z > -u^h/ \og{l/h). 

In this and the following sections we will assume that z = 0{h) so that we can work at 
a fixed energy level. That means that 

(4.3) P{z) = P-zQ + OHi^L<h^), 

P and Q are self-adjoint, and where Q = q'^{x, hD) G \I'^'°(X) is elliptic and has a positive 
symbol in a neighbourhood of T^j^X fl p~^([— 5, 5]). The estimates are uniform when we 
shift the energy level within | Re 2;| < 5q and hence we obtain the estimates in Theorem [T] 

For simplicity of the presentation we assume that T± have global defining functions, that 
is that F-i- are orientable. The only object that needs to be globally defined, however, is 



the escape function G given in (4.6). That involves only squares of defining functions, that 



is the (i(«, F-t)^, near and these are well defined and smooth. 
We start with the following 

Lemma 4.1. Let ip± he any defining functions ofT±: 

r± = {p : </?±(p) = 0}, rf(/?±fr±^0. 

Then, there exist c± G C°°(T*X;]R) such that 

(4.4) Hp(p± = =Fc|v?± , c± > zn neigh(i^'o), 
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and we can choose the sign of (f± so that 

(4.5) {ifi+,ip.}\K> co> 0. 

Proof. Since Hp is tangent to T± we have Hpip± = a±ip± and Hpip\ = 2a±ip\. To see that 
=Fa± > 0, we need to check that 

d 



Hpd{;T^) 



dt 



exp{tHpyd{;T^) 



\t=o 



in neigh(iro) 



But this follows from the assumption (1.5) which implies that 

d{exp{±tHp){p),T±f < Cexp{-9t)d{p,T±f , 0<t<T, 

for p in a T-dependent neighbourhood of Kq - see [35l Lemma 5.2]. 

To see (4.5) we note that dip±{p), p G Kq, are linearly independent and vanish on 
TpKQ C Tp(T*X) which is a symplectic manifold of codimension 2. Hence {H^^)p are 
linearly independent and transversal to TpKo, and 

because of the non-degeneracy of u, the symplectic form. If necessary switching the sign 
of one of the ip± we can then obtain (4.5). □ 

We define 

^l{p) + h/h 



(4.6) 



G{p) = xip) H 



^Hp) + h/h 



^ + Ci lo^ 



Xi{p)Gi{p) 



where: x ^ C^{T*X) is supported near Kq, with x = 1 on the set V in (3.25); Gi is 
described in Q xi e C^{T*X), 

Xi(p) = 1 , pe p-\[-26, 26]) n T*(o,2^)X ; 

supp Vx C {xi = 1}; and Ci is a large constant. Writing G"^ = G"'{x, hD) we observe that 

(4.7) WG'^uWjjk <\og{l/h)\\u\\L2 , \/k. 

We also recall an elliptic estimate: 

1 



(4. 



[P - iW)u\\l^ > 11(1 - ^r)(^ - ^W)u\\ > ^11(1 - i^2)n\\Hl - 0{h^)\\u\\L^ 



where Tpj G C^{T*X] [0, 1]) are as in Proposition 3.1 In fact, if Ei has the properties given 
in that proposition, 

11(1 - i^rMm = 11(1 - i^r)Ei{P - iW)u\\Hf + 0{h^ 



\u\ 



;i - i^^)E,{l - i,^){P - iW)u\\hi + ^(^' 



\u\ 



< C\\il - ij1'){P - iW)u\\l2 + 0(/i^ 



u 



which is (4.8). 



RESOLVENT ESTIMATES FOR NORMALLY HYPERBOLIC TRAPPED SETS 25 

The elliptic estimate shows that we only need to prove 

\\{P{z) - iW)u\\ > ^ 



log(l//i) ' 
for u satisfying 

(4.9) x'"{x,hD)u = u + OHk{h'^), \\u\\ = l 



where x has properties of, say, in (4.8 ). That is because the commutator terms appearing 
after this localization can be estimated using (4.8). 



Hence from now on we assume that u satisfies (4.9) with the support of x in a small 
neighbourhood of the energy surface p~^(0). 

We now proceed with the positive commutator estimate. Let Mq > 0, 3 z = 0{h), 
and calculate 

- 2 Im((P(^) - iW)u, (G"" + Mo \og{l/h))u) 

= -2lm{{P - zQ - iW)u,{G'" + Molog{l/ h))u) - 0{h^)\\uf 

(4.10) = ~i{[P,G'"]u,u) + 2Molog{l/h){Wu,u) + 2{Wu,G'"u) - 0{h^)\\u\\^ 

> h{{HpGYu, u) + 2Mo \og{l/h){Wu, u) - 2\\Wu\\ - 0(/ii/ii)||?if 

> h{{HpG)'"u,u) + Molog{l/h)\\w'^u\\^ - 0(/it/it)||Mf , 

where we used the fact that <W < y/W and chose Mq large enough. 

To analyze {HpG)'^{x, hD) we proceed locally using the invariance properties described 
in §3.3[ the resulting errors are of lower order. To keep the notation simple we write the 
argument as if ip± were defined globally (which is the case when r± are orientable). 

The crucial calculation is based on Lemma 14.11 

^f^ = -T^^ + -T^^ + ^o + ^l^^^ inneigh(Ko); 

ipj^ + h/h ip_ + h/h ^ 

here Rq is the term arising from Hp{x) and Ri from Hp{xi)- 
Put 



+ h/h 



We now recall the properties of Gi enumerated in §3.6t note further that supp Rq C 
{HpGi > 1}, hence for Ci ^ we may absorb the Rq term into the term arising from 
HpGi, and obtain the following global description of HpG : 

(4.11) HpG = ^1 + ^1+ Ri + Gi\og{l/h)a, 

where a G S(T*X), and 

a(p)> 1/2, d{p,K)>e>0, p G neigh(p-^(0)) , p G f/a . 
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We should now remember that using the rescahng (3.6) we are now in the semiclassical 
calculus with the h Planck constant. That means that the Weyl quantization is equivalent 
to the h quantization. 



Then (4.11) and the fact that we are using the Weyl quantization show that 

1 



{HpGY{x, hD) = $^ + $1 + Ci \og{l/ h)a'"{x, hD) + R^ + 0^^ (h' 
We now write 



so that, without writing the terms involving a'^{x, hD) and Ri, 
{{HpG)'"{x,hD)u,u) > {{<I>1 + <l>l)u,u) - OCh^)\\uf 

> ||$«||^+(^[<|.+,$_]n,n)-0(/.2)||«f 

(4.12) 

> Mh\\<!>u\\l + h{{(p+,lp.}'"{x,hD)u,u) -Oih^)\\uf 

> {{Mh{^l + ^^_) + h{^+,^.})'"{x,hD)u,u) - OCh^)\\uf 
where M is some large constant. Putting 



we calculate 



(1 + ^i)2(l + <^^)2 (1 + ^^) 2(1 + (^2)2 



(1 + ^1)^1 + )i (1 + + <^2_)i 

Hence 

satisfies ^ 

(p G hSi , 



and, using (4.5), we obtain near Ki^ 



0; 



[l + (P+)^{1 + (fit)-^ 2 



> h I Mi^l + f_) + r° - iih/h)r 

(1 + (y9^)2(l + (^1)2 2 

> ci/i , ci > . 
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We now return to (4.10) which combined with (4.7),(4.11 ), and the above definition of 
(p gives, for some large constant Mi, and u satisfying ( 4.9[ ), 



Ml log(l//i) II {P{z) - iW)u\\ \\u\\ > ((/i^"' + hRI' + Ci \og{l/h)a'" + Mq logil/h)W)u, u) 

> h{{(p'" + + \og{l/h)V")u, u) 

where, as > 0, 

h '^^ Cia + MqW > ^ r (x, hD) > -Ch , 
with the imphcation due to the sharp Garding inequahty. We also observe that 

hSi 3 (p + hh> cih, ci > , 



near p ^{{—6,6)). Furthermore, since u is assumed to satisfy (4.9), and as we have 
0{h°°) on such distributions, we obtain 

Milog{l /h)\\{P{z) - iW)u\\\\u\\ > h{{^''" + hb'")u,u) - 0{hHog{l/h))\\uf 

> C3hh\\u\\'^ , C3 > , 



which proves the bound (1.1) for Im^; = 



4.3. Estimates for Im^; > —v^h. To prove the estimates deeper in the complex plane we 
will use exponentially weighted estimates which use the same escape function G given in 



(4.6). We start with a lemma which is based on |39l Proposition 7.4]: 



Lemma 4.2. Let G he given by (4.6) above. Then for p, p' in any compact neighbourhood 
of Kq we have 

N 

<C( ) , iV>0. 



expG'(p) ^ I P- p' 



exp G{p') 



\{h/h)\ 



In particular, 



m{p) =^ exp G{p) 



is an order function for the \E'i calculus, that is, satisfies (3.20). 



Proof. For the reader's convenience we recall the slightly modified argument. We first claim 
that 



(4.13) 



'^±{p? + h/h ^ ^ I p-p' 
^±{p'f + h/h- '\{h/h)l 
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Since ~ d{», T±y, we have 

<^±(p)' + h/h < c(rf(p, r±)2 + h/h) < cw, r±)2 + \p' - p|2 + h/h) 

<C\v±{p'? + h/h+\p'-p?) 

= C'iMpT + h/h + {h/h){{p - p)/ih/W) 
<2C\ip^{pr + h/h){{p-p')/{h/hff. 



which proves (4.13). In other words, for 



G{p) = lo, 



dcf, ip'i{p) + e'^ 



h\ 2 
h. 



we have 



For x^C, 



Also, 



oo 
c ' 



\G{p)-G{p')\<C + 2\og{{p~p')/e). 

\x{p)G{p) - x{p')G{p')\ <G\p- p'\ log(l/6) + Clog((p - p')/e) 
\Xi{p)Gi{p) - Xi{p')Giip')\ <G\p- p'\ log(l/e) , 



with Gi as in ^3.6 thus to prove the lemma we need 

|p - p'l log ^ < C log((p - p')/e) +G, p,p'eQ(^ 

If we put if; = \p — p'\/ (Ge), this becomes 

, 1 log(t) + 1 „ 1 
e t e 

which is acceptable as the function t (log(t) + l)/t is decreasing. 



o2n 



□ 



We now consider {P{z) — iW)sG defined by (3.22) using this weight function G. Then 
using (3.24) and Lemma 3.2 (to understand adsc™ P), 

P[z)sG = P- ish{HpG)'"{x, hD) -zQ + O^^ {s^hh + sh'^U + h^) . 

2 

and 

= W + tsh\og{l/h){Hp,G,Wr{x,hD) + 0^^{s\hlog{l/h)y), 
where Gi and pi are as in (4.6). Hence, 

- lm{{P{z) - tW)sGU, u) = {{HpG)"" + W -Imz Q)u, u) 



For u satisfying (4.9), s > small, Im^ > —uoh for a sufficiently small uq, we can now 
proceed as at the end of to obtain invertibility: 

Cihh\\u\\ < \\{P{z) — iW)sGu\\ , Im z>—CQhh, \z\<Gh. 
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Since 

exp(±sG"'(x, hD)) = OL^^L'z{h-^) , 

that means that 

h'''\\u\\<Ci\\{P{z)-iW)u\\, lmz>-cohh, \z\<Ch. 



4.4. A global estimate. Here we show how the assumption (1.4) part (ii) oi Proposition 
3.1 give a global estimate; recall that the estimates of ^4.1 4. 3| applied to u supported in 
Xq. We fix a partition of unity on the interior of Xq 

1 = Xo + X? 

such that Xo = 1 on U2, supp Xi C {W = 1}, and with supp Xi C {W > 0} for i = 1,2. 
The results of ^ 4.1||4.2 and 4^ show that, in the notation of §1.1[ 

Im z , Im 2; > , 

(4.14) j{z,h)\\xou\\<C\\{Po{z)-zW)xou\\, ^{z, h) = \ h/ logil/h) , lmz = 0, 



h'' , Imz > — z/o/i . 



and, since XiW^ = 1, 

(4.15) collxi^^ll < \\{P{z)-tW)xiu\\, 



as implied by the hypothesis (1.4). 
Now, writing P{z) = P{z) — iW, 

\\p{z)ur = \\xoP{z)ur+\\ximur 

> \\P{z)xour + \\P{z)xiur - 11 [Xo, P{z)]uf - II [xi, P{z)]u\ 

- 2 {\\xoP{zM\\[xo,Pi^)M + \\xiPiz)u\\\\[xuPi^)]u\ 

> ||P(^);^o«f +||i5(^);^^^||2_2C'(||[;^,,p(^)]^f +||[^^,^ 

- wmur/c 

Since on the support of the commutator terms W = 1 and P{z) = Po{z), we have obtained 
Co\mz)-^W)ur > UPoiz) - ^W)xour + || (P(^) - ^)xi«f 

- Cidl [xo, (Poiz) - ^)]ur + II [Xi, (Poiz) - Ol^lH . 



Using (ii) of Proposition 3Jlwe obtain 



II [x„ (Piz) - ^)]ur < Chm^iPoiz) - ^)ur - Oih^)\\u^^' 



2 



< ChmP{z) - tW)uf - 0{h°^)\\u\\l , 

where ip ^ ^^(-^o) satisfies 

fsupp fsupp dxj — 1 • 
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We apply this estimate, (4.14), and (4.15), to get 

c^imz) - zw)ur > uPoiz) - iW)xour + \mz) - ^)xlur - o{h^)\\ur 

>^{z,h)\\xou\\' + co\\xiur-0{h'-)\\uf 
>^{z,h)i\\xour+\\xiuf)~0{h^)\\ur 
>{^{z,h)/2)\\ur. 
which completes the proof of Theorem [T] 

5. Results for resonances 

Here we briefly indicate how the proof presented in ^ adapts to give a resonance free 
strip. First we need to make additional assumptions on the operator guaranteeing mero- 
morphic continuation of the resolvent. 

Suppose that X is given by (3.1 ) with > 1. For simplicity we will assume that = 1, 
with obvious modifications required when for > 1. 

We make the same assumptions^ as in [39l (1.5)-(1.6)] and [3ll §3.2]: P = P{h) = P{h)* , 
P{h) = p'"{x, hD) + hp^{x, hD- h), Pie S^'\T*X) , 

(5.1) \^\>C ^ pix,0>{OVC, p = E^dp^O, 
3 i?o, V M G C°°(X \ Xo) , P{h)u{x) = P^{h)u{x) , 

where in X \ Xq = M" \ fi(0, R) 

(5.2) Poo{h)= Y,cio.{x-h){hDS , 

\a\<2 

with ) independent of h for | h) G C^i^"^) uniformly bounded 

with respect to h (here C^{W^) denotes the space of C°° functions with bounded derivatives 
of all orders), and 

^ aa(a;)r > (Vc)|^|^ "i^eW , for some constant c> 0, 

l"l=2 

(5.3) 

aa{x] h)^"" — y , as |x| — )■ oo, uniformly with respect to h. 

|q|<2 

We further take the dilation analyticity assumption to hold in a neighbourhood of infinity: 
there exist 9o G [0,7r), e > such that the coefficients aa{x; h) of Pooih) extend holomor- 
phically in x to 

[ruj : w G C" , dist(a;, S") < e , r G C , \r\> Rq, argr G [-e, + e)} , 
with ( 5.3[ ) valid also in this larger set of x's. 



assume that pi is of order 1 in ^ to make the case oi h — 1 easier to state. 
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We note that more general assumptions are possible. We could assume that X is a 
scattering manifold which is analytic near infinity and satisfies the conditions introduced 
in 



Theorem 2. Suppose P is an operator satisfying the dilation analyticity assumptions above 
and such that P{z) = P — z satisfies the assumptions of Theorem [7| Then for any x ^ 
C'^{X), x{P — z)~^x> continues analytically from Imz > to Imz > —Voh, \z\ < 5o, and 



C^h~Hog{l/h) , lmz = 0, 



(5.4) \\x{P - z)-W\l2^l^ < 

for \z\ < 6q. In other words, there are no resonances in a strip of width proportional to h. 



Sketch of the proof: The proof follows the same strategy as the proof of the estimate 
for Im^; > —v^h in Theorem^ but with W replaced by complex scaling with angle 



6 h log(l / h) . That requires a finer version of Lemma 4.2 which is given in [311 Proposition 
7.4]. In particular, the choice of the cut-off function xi has to be coordinated with complex 
scaling (see also [391 §4-2]). The same exponential weight can then be used, following 
the arguments of [391 §8-4], but without the complications due to second microlocalization 
needed there. 

This provides the bound 0{h~^) for the norm of the analytically continued cut-off re- 
solvent, x{P ~ for > —vqH. To obtain the bound on the real axis we can 
proceed either as in §4.21 using the "semiclassical maximum principle" - see for instance 
[SI Lemma 4.7] or [TUl Lemma A. 2]. □ 

Ideas used in the semi-classical case provide results in the case of the classical wave 
equation. We first note that if P = P(l) satisfies the assumptions above then the resonances 
are defined as poles of the meromorphic continuation of {P — \^)~^ from Im A > to Im A > 
— co|ReA| - see [36]. When Poo = —A and the dimension, n, is odd, the meromorphic 
continuation extends to the entire complex plane (that is why we use the parametrization 
z = A^, and when n is even we pass to the infinitely sheeted logarithmic plane) - see [37] . 
Theorem |2] implies that for x ^ C^(^)) 

(5.5) ||x(-P - A^)"^xlU2^L2 < C^|A|'' , ImA > -ai , I Re A| > ao , a,/3>0. 

To relate this to energy decay we procceed in the spirit of |9j. Suppose that the operator 
P satisfies the assumptions above with h = 1 and consider the wave equation for P with 
compactly supported initial data: 

(5.6) {D^ - P)u = , u\t=o= uq , Dtu\t=Q= ui , supp Uj dV X . 



The local energy decay results are different depending on finer assumptions on P which we 
state as three cases: 
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Case 1 


P 


K"\_B(0,i?.()) = 


R"\B{0,Ro) 


n odd 


Case 2 


P 


M"\B(0,flo) = ~^ 


R"\B{0,Ro) 


n even 


Case 3 


P 


M"\B(0,Ro) — Poo\ 


M"\B(0,i?o) 


any n 



where Poo is an elliptic operator close to the Laplacian at infinity - see (5.2) and (5.3) - 
with h = 1. 

Theorem 3. Let P be an operator satisfying the assumptions above with h = 1. Let 
U,V G X be bounded open sets, and let \E' G C°° (M) be an even function such that 



(5.7) ^(a;) = l 



for X G M m cases 1 and 2 



(x) = near in case 3. 



for X > 1 in case 3 

Suppose that P has neither discrete spectrum nor a resonance at 0. Then there exists K > 



such that the solutions of (5.6) with 

ll'^^oll/fA'+i < 1, ll'^^lllz/A' < 1, '^{VP)Uj = Uj 



satisfy the following local energy decay estimates: 
(5.8) 

where the constant C (Cm) depends on U and V (and M) only 



Cexp(— at) , in case 1, 

{\u{t,x)\^ + \dtu{t,x)\'^)dx < { Ct-"+Mogt, incase 2, 

CMt'^^ , VM > 0, zn case 3, 



Proof. We first note that it is enough to obtain the estimates xU{t)x '■ — > L^ where 
X G C^{X) and 

^^^^de^sintVP 



p 

To do that we follow the standard procedure (see |ll],[9l §4] and reference given there) and 
perform a contour deformation in the integral: 

• r+oo 

(5.9) xu{t){p + t)-^/'^{Vp)x = / e-^''xim - Ri-m^' + tr^'/'^wxdx , 
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for t > 0. The contribution of R{—X) in the spectral projection can be ehminated by 
contour deformation when t > - see [Hi Sect. 4]. Hence 

r>+00 



(5.10) xump+^r''^''^i^)x 



2^ 



-it\ 



xi?(A)(A2 + i)-^/2^(A)xc^A, t>0. 



In case 1, i.e., odd dimensions and P = —A in the exterior of a (large) ball, we use the 



27, < 7 < «!. This gives (5.8) in that 



estimate (5.5) to deform the contour to F = I 
case. 

In the case of a compactly supported perturbation of —A and n even, we have to modify 
this argument because the resolvent has a branching point at A = 0. Thus we deform the 
contour near to 

{X = X — icix, X > 0} U {z = X + icix, x < 0}, for ci > 0, small. 

We use the usual estimate for the resolvent near 0: 

\\xRiX)x\\<CM\X\"-'\iogX\ 

in any sector | argA| < M - see for instance ^S', §3]. The dominant part of the integral 
(5.10) comes from the contour near which gives 



„n-2 ■ 



X 



logxe-'^'dx < Ct""+Mogt 



which is the estimate in case 2. 

For case 3, that is the case of ^ 1, we consider the analytic extension of that function, 
^, with the property that d"^ = 0{\ lmz\°°) (the defining property of the almost analytic 
extension - see [16, Chapter 8]) is supported in a set where P has no resonances - see FigjT| 
We deform ( |5.10[ ) to a contour which for \z\ > 1 is the same as before, and for \z\ < 1 is 
as in Figjl] By Stokes's formula we get exactly the same contributions as in case 1 (since 





^ = 1 




support of 



Figure 1 . The contour deformation in case 3 and the support properties of 
the almost analytic extension of ^P. 

near 0, = 0) with an additional term 



(5.11) 



27r 



(z) e-'''xiR{z){z' + iy'^x dz 
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where Q is the support of between the real axis and the new contour (shaded in Fig. 
Since (z) = O {\ lmz\°°), a repeated integration by parts shows that this last term is 
O {t~°°) (in the energy norm). □ 



Proof of Corollary^- We follow the argument of Burq [S]. The left hand side of (5.8) is 
bounded by the same quantity at t = 0, and in particular by ||mo||^i + ||wi||^2- The estimate 
(5.8) shows that, in case 1 (that is, the case considered in Corollary]^ 
by e~'^*||Mo||_H-if+i + Interpolation between these two estimates gives (1.2) 



it is also bounded 

□ 
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